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nomenclature 


B A constant in Eqs. (2) and (4) 

C A constant in Eqs. (2) and (4) 

C Negative of C 

K Thermal conductivity (WATT/cm-°k) 

P Phase volume fraction 

Pj Porosity 

P Mass fraction 

m 

R Thermal resistance 

T Absolute temperature (°k) 

E Modulus of Elasticity of solid particle 

g Gravitational constant 

z Depth below the surface 

p Density 

p - Average density 

u Poisson's ratio 

Subscripts 

c Continuous phase 

d Discontinuous phase 

eff Effective 

S Solid 
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INTRODUCTION 


A comet is generally regarded to be composed of three primary regions: 
the nucleus or Kemal, the coma (a plasma region surrounding the nucleus), 
and the tail. Since not much is known about the composition and structure 
of the nucleus, it is thought that some useful information may be inferred 
by investigating the heat and mass transfer of the nucleus. Before con- 
ducting these investigations, it is necessary to study the material thermal 
properties of possible nucleus models. 

The object of the present study program is to perform the following: 
Part I: To recommend a value for the thermal conductivity of the frost 

layer and a value for the thermal conductivity of the water-ice 
solid debris mixture. The basis for the recommendation is a 
nucleus model consisting of water-ice mixed with basaltic and 
meteoritic materials in the form of dust and agglomerated par- 
ticulates, all of which is surrounded by a layer of water- frost, 
assuming a constant thermal conductivity in each of the layers. 

Part II: To recommend a value for the thermal conductivity of the porous 
structure of Part II as a function of depth only. The model for 
this part is the same model as in Part I with the following addi- 
tional assumption. The assumption is that the frost layer and 
water-ice have been sublimating and evaporating, leaving a porous 
structure composed of the solid basaltic and meteoritic material 
with residual gases partially confined in the porous structure to 
varying amounts over some depth. This conductivity function should 
be regarded as a simplified approximation to be used for initial 
calculational estimates of heat transfer in the nucleus. 
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COMETS CONSIDERED 

i) Commet Encke: With a minimum temperature of 140°K, Radius: 4KM. 

ii) Comet Halley: Min, Temp: 49?K, Radius: 15KM. 

ANALYSIS FOR PART I 

For the Water Ice-Solid Debris Mixture, the Cheng and Vachon Model [1] 

has been used to calculate the approximate effective thermal conductivity 

of this layer. The nucleus of a comet normally consists of 70% or more 

♦ 

of ices by mass [2] and the rest, basaltic and meteoritic material. As a 
first approximation, in the absence of detailed information on the composition 
of the meteoritic material (other than that it might contain Fe, Ca, Mn, 

Mg, Cr, Si, Ni, Al, etc), the mixture of basaltic and meteoritic material 
has been approximated first with the properties of basalt and then with 
properties of basalt and iron which are mixed in different proportions. 

(A) : When the solid debris has been approximated with the properties 

of basalt, the above model is simplified to a mixture of water ice and 
basaltic material. Since the temperature of the nucleus varies with its 
distance from the Sun (Example: At Aphelion [3] temperature of nucleus 

of Encke = 142° K and that of Halley = 49°K),the thermal conductivities 
have been calculated for the various assumed temperatures of the nucleus. 

The calculations have also been performed with varying proportions of 
ice and basaltic material. 

(B) : When the solid debris was approximated by a mixture of basaltic 

material and iron, the model is simplified to a mixture of ice, basaltic 
material and iron. Again the calculations have been performed for different 
assumed temperatures of the nucleus (variation of thermal properties within 
the nucleus for a given average temperature have not been considered) and 

for varying proportions of water ice and solid debris. The calculations 
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have also been performed with varying proportions of basalt and iron for 
a given ice, solid debris ratio. 

A linear variation of thermal conductivity with temperature has been 
assumed for the water ice and using the values given in [4], The follow- 
ing equations were obtained for the thermal conductivity and the desnsity 
variation of ice with temperature 


a) K. ce = [0.0481156 - 9.2444 x 10' 5 T]Watt/cm-K 

b ) p ice = [0.94323 - 9.2444 x 10" 5 T] gm/cm 3 T in °K. 


( 1 ) 


The following property values were 
Basalt : p * 2830 Kg/m 3 

E = 2.2 x 10” N/m 2 

Iron: p = 7890 Kg/m 3 

E = 2.07 x 10” N/M 2 


used for the Basalt [5] and Iron [6] . 
K s = 1.34 watt/m°K 
v = 0.2 

K s =63.7 watt/m°K. 
v = 0.3 


K. is of the order 3.8 watt/m°K at 100°K 

ice 


1(A) 

For this model, we see from the above that for the temperature 

K. 

1 C 6 

range of 30-200°K, the ratio y is in the range 1 to 4. Hence, 

solid 

from [1], for the case K^ > we have the effective or equivalent thermal 
resitance given by 


p = 2 tan" 1 B /r 

eff /(- c’ (k d - kj IK. + B (K d - K c )] ) 2 Vi 


-C' (K. 


K ) 
c 


K + B(K. - K ) 
c d c J 


} * 


1-B 


( 2 ) 
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Where K = Thermal conductivity of continuous phase 

c (water-ice in this case) 

K, = Thermal conductivity of discontinuous phase 
(basalt in this case) 

= Discontinuous phase volume fraction 


B = /Ip 
C = -4/B 
C’ = -C 


d/2 


Then the approximate effective thermal conductivity is given by 
K 


eff 


1 _ 

R 


eff 


(3) 


Thus, the has been calculated for different temperatures of 

nucleus and for different phase volume fractions of the discontinuous 
phase (basalt in this model) . 

ICB) 

In this model, we have a three phase mixture of water-ice, basalt 
and iron and we have to find the approximate thermal conductivity (K^^.) 
of this three-phase mixture. Cheng and Vachon [1] give the following 
method for calculating 

Since the major portion of the comet nucleus is made up of water- 

ice, we take water-ice to be the continuous phase in applying the equation 

developed by Cheng and Vachon for 3 phase mixtures. 

Let K = Thermal conductivity of the continuous phase (i.e. 

L of water ice) 

K. = Thermal conductivity of the first discontinuous 
d l phase (i.e. of Fe) 

K, = Thermal conductivity of the second discontinuous 
°2 phase (that is of basalt) 
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Thus, from the given range of values for the above, we have 


K 


d 


1 



> 10 



< 1 


In order to determine the effective thermal conductivity (K^^) 
of this 3 phase mixture, the 3 phase mixture is considered to be 
reduced to a two-phase mixture in which the continuous phase is composed 
of the original continuous phase of the three phase mixture and the 
second discontinuous phase. The discontinuous phase of this two phase 
mixture is the first discontinuous phase of the three-phase mixture. 

For this two-phase mixture. 


P - P 

<1 ci. 


P = p + p , 
ce c d. 


K d " \ 

K = Thermal conductivity of the effective continuous phase 
00 

P = Phase volume fraction of the effective continuous phase 
ce 

K ce , the effective thermal conductivity of the original continuous 
phase mixed with the second discontinuous phase of the 3 phase mixture, 

is determined as per the method given in Part I (A) . 

K, 


Now, since 


K 


> 10, from ref [1], the effective resistance R 


ce 


eff 


is given by 


Wf TIcmTT T~k ) [K + B' (K, He HT 
cL ce J 1 ce I d, ce^ J 


In /{^ce 


„ ,K, K A B. /{C ' (K , - K )> 

+ B 1 ( d l - ce) j + 1 d 1 ce J 


/<K ce * B 1 <*<1, - K ce» - B l - Kj) * 1 - B ! 

1 2 1 K 

ce 


( 4 ) 
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where B, = /3P, /2 
1 d l 

C» = -C = + 4/Bj 

and the effective thermal conductivity is then given by 

K eff = (5) 

K eff 

This calculation has been performed for different assumed temper- 
atures of the nucleus (Here also, the thermal property variation inside 
the nucleus has not been considered) and for different phase volume 
fractions of the continuous phase. For a given value of phase volume 
fraction of water ice, calculations have been performed for varying 
proportions of basalt and iron. 

ANALYSIS FOR PART II 

For the model in Part I, we now assume that the frost layer and 
the water ice has been sublimating and evaporating which occurs as the 
comet nears the sun (distance less than 2 AU) , thus leaving a porous 
structure . 

A simplified approximation for the effective thermal conductivity 
as a function of depth is obtained using the theoretical model developed 
by Khader and Vachon [5] for heterogeneous mixtures. In the application 
of the above model to our case, it has been assumed that since the gas 
is at low pressure the thermal conductivity of the residual gases is 
quite small compared to the thermal conductivity of basalt and meteoritic 
material. Hence, the thermal conductivity through the void can be 
neglected. Also, since the temperatures involved are of t-he order of 
100°K, radiation through the void space (i.e. through the gases at very 
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low pressures) could be neglected as being very small in comparison to 
conduction through the solid. 


11(A) 

Thus, with these simplifications, from Ref [5] we get the following 
expression for the effective thermal conductivity when the solid is 
approximated with the properties of basalt. 


K eff = 3(1 - P,)(l - cj Z 2/3 ) K s 


[J ( 3 1/ 3 - 1) + (4 - 1.2 C 2 Z 2/3)] 

2 c ^ Zi 


( 6 ) 


_ ? ■ L/ J 

where C ( = 2 g ptr(l-v ) 

16 H (1 - P,) 

p = Average Density of the Solid. 

P, = Porosity of the solid material. 

The calculations have been performed for different porosities and 
for different densities up to a depth of 2500 meters. 

11(B) 

When the solid material is approximated by a heterogeneous mixture 
of basalt and Fe, the problem is solved in 2 stages. First, an effective 
value of thermal conductivity is obtained for the heterogeneous mixture 
of basalt and Fe using the Cheng and Vachon model [l]. 

Then the model for the second stage consists of the solid material 
whose thermal conductivity was determined in Stage 1 and whose voids are 
assumed to be filled with low pressure residual gases. The effective 
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thermal conductivity of this stage is determined using the Khader and 
Vachon model as explained for 11(A). 
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RESULTS AND DISCUSSION: PART I 

Part I (A) : 

Fig. 2 shows the variation of effective thermal conductivity with the 
average temperature when the solid debris is approximated with the 
properties of basalt for different values of the density of the dis- 
continuous phase at a given value of mass fraction of the discontinuous 
phase [For example: The above graph is for the case when there is 80% 

of water-ice and 20% of basalt by mass] . This figure shows that the 
effective thermal conductivity of this layer is of the same order as 
that of water-ice and that its variation with temperature is also linear. 
The reason for this is that the major portion of the nucleus (about 70-80%) 
is made up of water-ice and the thermal conductivity of basalt is also 
of the same order of magnitude as that of water-ice. The figure also shows 
that at a given temperature there is not much variation in with change 

in density for this case, especially near temperatures of the order of 200°K. 

Fig. 3 shows the variation of with mass fraction of solid debris 

(i.e., discontinuous phase approximated by basalt in this case). It shows 
that the variation is not appreciable in the range of mass fraction pre- 
dicted for the nucleus of a comet and that there is a slight decrease in 
the thermal conductivity as mass fraction of solid debris is increased. 

It also shows that the variation is more pronounced at temperatures near 
40°K than near 200°K where it is almost flat. 

Part 1(B) : 

For the model of 1(B), variation of thermal conductivity with the 
total mass fraction of basalt and iron is shown in Fig. 4 for different 
proportions of Fe and basalt. As can be seen from this figure, though 
the curves start at almost the same value of for very low values of 

mass fraction (the reason for which is that practically the whole nucleus 
in this case will be made up of water-ice), the slope of the curve changes 
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from being negative for the model with only basalt, to positive for the 
model with only Fe as the solid debris material. The reason for the 
positive slope in the latter case is that since the thermal conductivity 
of iron is much greater than that of water-ice, the thermal conductivity 
of the mixture increases as the percentage of Fe increases. When Fe and 
basalt are in the ratio 1:1, even though the slope is positive, the 
increase is not marked. 
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PART II 

Figs. 5, 6, 7 show the plots of variation of the effective thermal 
conductivity with depth for different values of porosity, each figure for 
a particular proportion of Fe and basalt. They all exhibit the same 
general trend of increase of with depth, the increase being sharper 

during the earlier part and then becoming gradual. The reason for this 
is that as the depth increases the contact area between the particles 
increases due to the increase in loading and consequently smaller thermal 
contact resistance and in turn thermal conductivity increases. As the 
porosity increases, the amount of solid material per unit volume decreases 
and hence thermal conductivity decreases with porosity. 

A comparative study of Figs. 5, 6, 7 is made in Fig. 8. Comparing 
curves 1 and 5 we see that the thermal conductivity increases by 
a factor of about 75 when basalt is replaced by Fe. This marked increase 
in with increase in percentage of Fe compared to basalt is because 

of the very high value of thermal conductivity for Fe compared to basalt. 
This figure also shows that the curve of sharply rises for the first 

few hundred meters and then becomes very gradual (almost flat) . Thus after 
about 500 meters, the effect of depth on thermal conductivity is not 
pronounced. This is because after this depth, no appreciable increase in 
contact area between the particles occurs with depth. 

Fig. 9 shows the variation of effective thermal conductivity with 
depth for the case 11(A) [i.e. with only Basalt as solid material] for 
the first 200 meters from the surface. This shows that even though the 
graph still follows the trend of Figs. 5, 6, 7, actually the major portion 
of the increase in the value of takes place within a few meters from 

the surface and afterwards the increase is very gradual. Since the same 
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pattern was followed for the other cases (for varying proportions of 
basalt and Fe) they have not been plotted. 

Fig. 10 illustrates the variation of with density at differ- 

ent depths when the solid material was approximated with the properties 
of basalt. The increase in thermal conductivity with density is very 
gradual and the reason for the increase in with density is the 

increased contact area between the particles which in turn is due to 
increase in loading due to heavier mass/unit volume at a given porosity. 
Fig. 11 also illustrates the same. 



Effective Thermal Conductivity K ^ (Watt/cxn-K xlO^j 
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FIGURE 2: EFFECTIVE THERMAL CONDUCTIVITY AS A FUNCTION OF TEMPERATURE 
AT CONSTANT DENSITY FOR PART 1(A).' 



Effective Thermal Conductivity . K (Watt/cm-^K xlO ) 


is 



P ra _ (Mass Fraction) 

FIGURE 3: EFFECTIVE THERMAL CONDUCTIVITY AS A FUNCTION OF MASS FRACTION OF 
SOLID DEBRIS (BASALT + METEORITIC MATERIAL ) AT CONSTANT TEMPER' TORE 
FOR PART 1(A). 
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p = Mass Fraction Of (Fe + Basalt) 
m 


FIGURE 4: EFFECTIVE THERMAL CONDUCTIVITY AS A FUNCTION OF MASS FRACTION OF SOLID 

DEBRIS (BASALT + METEORITIC MATERIAL ) AT CONSTANT TEMPERATURE FOR PART I (B) . 






Effective Thermal Conductivity K (Watt/cm-?K x!0 J ) 
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Depth Z in meters 


FIGURE 6: EFFECTIVE THERMAL CONDUCTIVITY AS’ A FUNCTION OF. -DEPTH AT CONSTANT POROS ITY 

FOR PART 11(B). 
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. IRON AS SOLID MATERIAL. 



Z in meters 


IOURE 7: EFFECTIVE THERMAL CONDUCTIVITY AS A FU! 
FOR PART 11(B). 




0 



Depth Z in meters 


FIGURE 8: COMPARISION OF EFFECTIVE THERMAL CONDUCTIVITY FOR DIFFERENT PROPORTIONS 

OF BASALT AND IRON AS A FUNCTION OF DEPTH AT CONSTANT POROSITY FOR PART 1 1(B) 



Effective Thermal Conductivity ^ e ff xlO ) 
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Depth Z in meters 


:"I r -URE 9: EFFECTIVE THERMAL CONDUCTIVITY AS A FUNCTION OF DEPTH AT CONSTANT POROSITY 

FOR PART 11(A). 
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FIGURE 10: EFFECTIVE THERMAL CONDUCTIVITY AS A FUNC1 
FOR PART 11(A). 
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Depth Z in meters 


.._i T) ' 

,, 

FIGURE 11 : EFFECTIVE THERMAL CONDUCTIVITY AS A FUNCTION OF DEPTH AT CONSTANT DENSITY 
FOR PART 11(A). 
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